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0. Introduction 

It has long been known for a class of exactly solvable models in statistical mechanics 
that the correlation functions are expressible in terms of Painleve transcendents or more gener- 
ally solutions to completely integrable differential equations possessing the Painleve property 
[WMTB] , [MTW] , [SMJ]. (For a review see [Mc].) It was the great insight of the Kyoto 
School to invent the notion of a r-function associated to a system of completely integrable 
differential equations [SMJ], [JMMS], [JMU], [JM] and to realize that in this class of 
solvable models the correlation functions are r-functions. 

In the theory of random matrices it has also long been known (see, e.g., [Po], [Mel]) 
that the level spacing distributions are expressible in terms of the Fredholm determinant of 
an integral operator. In a particular scaling limit the kernel of this integral operator is the 
sine kernel. For this case it was shown in [JMMS] that the Fredholm determinant was also a 
r-function. This then leads in the simplest case to formulas for the level spacing distributions 
in terms of a particular Painleve V transcendent. More generally, the deformation equations 
consist of a compatible set of nonautonomous Hamiltonian equations. 

Since the original work of [JMMS] on random matrices and impenetrable bosons, the 
methods for deducing the relevant deformation equations and properties of the r-functions 
have been simplified and generalized, both in the study of random matrices [BTW] , [Me2] , 
[Dy], [MM], [TW1], [TW2], [W] and quantum correlation functions [IIKS], [IIKV]. The 
Hamiltonian structure of the underlying deformation equations can most simply be understood 
as a nonautonomous version of the equations induced by spectral invariants on a loop algebra 
with respect to a suitably chosen classical R-matrix Poisson bracket structure [H]. 

In the following section, the problem of computing level spacing distributions for the 
Gaussian Unitary Ensemble (GUE), both in the "bulk of the spectrum" and at the "edge 
of the spectrum" is formulated in these terms — the two cases differing only in the kernels 
resulting from the different scaling limits. For the first case, this is the above mentioned sine 
kernel 

1 sin7r(x — y) 
tv x — y 

while for the second it is the Airy kernel 

AiMAi'frJ-AifrJAi'fr) 
x-y 

In §2, the resulting Hamiltonian equations are explained in terms of spectral invariants on 
the dual space si (2) of the loop algebra si (2), with respect to a standard classical .R-matrix 
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Poisson bracket structure, restricted to Poisson subspaces consisting of rational coadjoint 
orbits. 

1. Spectral Distributions for Random Matrices 

la. Fredholm Determinant Representation 

The GUE measure on the spectral space of N x N hermitian matrices with eigenvalues 
(X 1 ,...,X N ) is [Mel] 

P N (X U . . . , X N ) dX 1 --- dX N = C N J[ \Xj - X k \ 2 exp I - J2 X j I dX i ■ ■ ■ dX N- (1-1) 

j<k V j ) 



Averaging over all but n of the eigenvalues gives the n-point correlation function [Mel] : 

A! f°° f°° 

Ru,n(Xi, . . . , X n ) = — — — — / • • • / P N (Xi, . . . , X N ) dX n+ i ■ ■ ■ dX 

yl\ — n). j-oo J -oo 



3,k=l 



where 

N-l 



' Vv 

(1.2) 

= det[K N (X J ,X k ) 



K N (X, Y) := J2 <Pk(X)<Pk(X) ■ (1-3) 

k=0 

Here {<pf-(X)} is the sequence obtained by orthonormalizing the sequence {X k exp(— A 2 /2)} 
over (— oo, oo). In particular R\^{X), the density of eigenvalues at X, equals Kn(X,X). 

For fixed X, 

Ri,n(X) ~ -V2N as N -> oo. (1.4) 

7T 

The scaling limit in the bulk of the spectrum at the point X is the limit 

A — > oo, Aj — > A , such that = (A^ — Ao)-Ri,at(Ao) is fixed. (1.5) 
In this limit the scaled n-point correlation functions R n (xi, . . . , x n ), defined by 

R n (x 1 ,...,x n )dx 1 ---dx n = lim R nN (X 1 , . . . , X n )dX 1 ■ ■ ■ dX n , (1.6) 

N — * oo , X j — »■ Xq 
xj fixed 

are given by [Mel] 



Rn(xi,...,x n ) = det[K(xj,Xk) (1.7) 
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where 

T w \ 1 sin7r(x — y) .„ n . 

K{x,y):= ^ y -l. 1.8 

TV x — y 

Now consider a set of "spectral windows" 

J = Ji U • • • U J m 

where 

{Jk = [«2fc-i, «2fc] C R}fc=i,... m , ai < a 2 < • • • < a 2m (1.9) 

is an ordered set of disjoint intervals on the real axis. The probability that a randomly 
chosen GUE matrix has precisely {ni, n 2 , . . . , n m } scaled eigenvalues in the respective intervals 
{Ji-, Ji-, ■ • • i Jm} is given by the level spacing distribution: 



E(n 1 ,...,n m ,J)= { ^ f,?^ (1.10) 
ni! • • -n m ! dA/ • • -dAm" 1 A 1= ...=A m =i 

with n = ni + • • • + n m . Here D(J; A) is the Fredholm determinant 

D(J;A) :=det(/-K) (1.11) 
of the integral operator K with kernel 

m 

Kj(x, y) := ^ K(x, y) X Jk (y), (1.12) 
fc=i 

where X T (y) is the characteristic function for the interval Jfc, and 

A:=(Ai,A 2 ,...,A m ) (1.13) 

denotes the set of generating function parameters. The case of one interval is in [Mel] . The 
more general formula (1.10) was derived in [TW1]. 

In the following, we shall see that D(J; A) is determined by a system of partial differential 
equations consisting of a compatible set of nonautonomous Hamiltonian equations defined on 
a suitable phase space. The quantity of central interest will be the total differential 



u{a) := d a \ogD{J- A) = £ dl °^ A ) d(Xi (1.14) 

i=i ° ai 

of the logarithm of the generating function D(J; A), taken with respect to the set of parameters 

a := (ai,...,a 2m ) (1-15) 
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giving the boundaries of the spectral domains. The dynamical equations of Jimbo, Miwa, 
Mori and Sato [JMMS], which determine the dependence of this quantity on the parameters 
(ai, . . . , a 2m ), are given below, following the approach developed in [TW1]. We note that in 
[JMMS], [TW1], and [TW2] the case of equal Aj's only was considered in the derivation of 
the equations. 

lb. Dynamics of Distribution Functions 

The key object in what follows is the resolvent operator 

R:= (l-K^K. (1.16) 

Introducing the functions 

Q(x,a) := (1 -K^A (1.17a) 

P(x, a) := (1 — K)~ l A\ (1.176) 

where 

A(x) := -sin7rx, (1-18) 

TV 

the kernel of the resolvent operator can be shown to be given by 

R{x,y)=y x k x jAy)i x rv, (1.19a) 

A — ' x — y k 

k=i y 

R(x, x) = J2 X k a ) p ( x > a ) " ^"fo a ) Q( x > a )) *j k ( x )- ( L1% ) 

fc=i ^ x x ' 

Using the identity 

d a log (det(J - £)) = -tr ((/ - ky x d a K} , (1.20) 

it is easy to show that 

2m 

u>(a) = J2(-l) j+1 R(*j,aj)da j , (1.21) 
i=i 

where R(aj,af-) is obtained by taking the limits (x,y) — > (aj,ak) inside J. Let 

qj := lim \/A7(2(:E,a) (1.22a) 

x—>aj 

Pi := Jim y/\~P(x,a), j = 1, . . . , 2m. (1.226) 



Then we have 



R( aj ,a k ) = qjPk m \ j^k, (1.23a) 
ctj — a k 



2m 



R(aj,a,j) = 7r 2 g| + pj +^2(-l) k R(a j ,ak)R(a k ,a j )(a j - a k ) ■ (1.236) 



fe=i 



Thus, the parameter dependence of det(I — K) is determined through the quantities 
{Qj,Pj}j=i,...2m- Differentiating equations (1.17a,b) with respect to the parameters cij and 
taking the appropriate limits, these may be shown to satisfy the following system of equations: 

dqj_ 
da k 
dpj_ 
da k 

H ^ „a-„, . , , 0-2 1) 

dpj_ 
daj 



(-l) k R(aj,a k )q k , j^k 
(-l) k R( aj ,a k ) Pk , j^k 



Pj ~ 5^( _1 ) fcjR ( a J' ak ^ qk 



This overdetermined system, which we refer to as the JMMS equations, is in fact Frobenius 
integrable, and has a Hamiltonian structure which we summarize in the next subsection. 

lc. Hamiltonian Structure of the JMMS Equations 

For convenience, define a new normalization of the coordinates 



i 1 
-x 2 j, <?2j+i := - 

P2j ■= -iy2j, P2j+i ■= y 2 j+i 



q2j := --x 2j , q2 J+ i := ~x 2j+ i ^ 



and let 

In terms of these quantities, we have 

2m 

u(a) = s ^Gj{x,y) daj. (1-27) 

3 = 1 

Defining canonical Poisson brackets: 

{xj,x k } = {yj, y k } = 0, {xj,y k } = S jk , h k = 1, . . . , 2m, (1.28) 
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the system of equations (1.24) is seen to have the nonautonomous Hamiltonian form 



d a Xj = {xj,u(a)} and d a y i = {yj, u(a)} . (1-29) 

Moreover, the set of Hamiltonians {Gj }j=i,...2m are easily verified to be generically functionally 
independent and in involution: 

{Gj,G k } = for all j,k=l,...,2m. (1.30) 

From this follows that the system (1.29), or equivalently (1.24), is Frobenius integrable, and 
that uj is locally an exact differential: 

U = <i a log T, (1-31) 

where the value of the tau- function r along the integral curves coincides with log (det(J-K)). 

The significance of this Hamiltonian structure in terms of spectral invariants on loop 
algebras will be explained in §2. First, another similar system occurring in the computation of 
level spacing distribution functions of scaled random matrices at the "edge of the spectrum" 
will be described. For the details of what follows see [TW2] . (Here also only the case of equal 
Aj's was considered, the generalization to the case of unequal A/s being straightforward.) 

Id. The Airy Kernel System and Distributions at the Edge of the Spectrum: 

Scaling at the "edge of the spectrum" [Mo], [Fo], [TW2], corresponds to choosing 
Xq ~ ±v / 2iV and gives rise to a different Fredholm determinant, in which the sine kernel (1.8) 
is replaced by the Airy kernel 



K(x,y) : - A (x)A'(y)-A'(x)A(y) 



= lim 



x-y 



(1.32) 



TV^oo 2 1 / 2 jV 1 /6" 

where now, A(x) is an Airy function 

A(x) = Ai(x). (1.33) 

The logarithmic differential of the Fredholm determinant (1.14) in this case is also given 
by formula (1.21) where, using the same notation as above, the resolvent kernel R(x,y) is 
still of the form (1.19a,b), with the functions Q(x, a), P(x, a) obtained by replacing A{x) in 
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eq. (1.17a,b), (1.18) by (1.33). Define {qj,Pj}j=i,...2m as in (1.22a,b), and introduce two 
further quantities: 



j 

Then similarly to the previous case, we have, 



u := ^ Aj / A(x)Q(x,a)dx (1.34a) 

j—\ Ja2j-i 
m „ a2j . 

v-J^Xj A(x)P(x,a)dx. (1.346) 

a— i Ja 2 j-i 



R(a^a k ) = qjPk m \ j^k (1.35a) 
a,j — a k 

R(aj, a 3 ) = ]T(-l) fc {QjP * ~_ P f k)2 + v) - a i?i 2 - 2 Pi <&« + 2gJ V , (1.356) 
Mi 



(ij — a k 



where all the limits are again taken within the intervals Jj, J k . Differentiating with respect to 
the parameters {aj}j=i,...2m gives the system 



dq 



da k cij - a k 

oa k aj - a k 

9qj ^r f uk q 3 Pk-p 3 qk 

— = - } (-1) J - — p k + ajqj +pjU-2qjV 

da j t^j <•:, 
du 



(1.36) 



dv 
daj 



{-l) 3 Pjqj. 



Again, introducing new coordinates {xj, Hj}j=o,...2m by: 

i 1 

?2j :— _ 2 X2j ' q2j+1 ''~ 2 X2j+1 

P2j ■■= -iytj, P2j+i ■= V2j+i (1-37) 

1 

u := y , v := -x . 

and defining canonical Poisson brackets 

{xj,Xk} = {Vj,Vk} = 0, {^-,2/*.} = ^, j, fc = 0, ...,2m, (1.38) 
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eqs. (1.36) have the Hamiltonian form 



d a Xj = {xj,u(a)} and d a yj = {yj,u(a)} , j = 0, ...,2m, (1.39) 

where 

2m 

uj(a) = d a \og(det(I - K)^ = ^Gj{x,y) da t (1.40) 

i=i 

and 

n< \ 2 1 2 , 1 2 1 ^ (xjyk - x k yj) 

Gj(x, y) = Vj - -ajXj - Xj-yj-yo + jXjX - - ^ — — • (1.41) 

fc=i j a fc 

It is straightforward to verify that these G/s are also in involution, implying the Frobenius 
integrability of the system (1.40) and the existence of a r-function as in eq. (1.31). However, in 
this case, the phase space is 2(m + l)-dimensional, so the G/s alone do not form a complete set 
of commuting integrals. It is easily verified, though, that the following additional independent 
function 

^ 2m 

G (x,y) :=y 2 -x --J2 x2 i ( L42 ) 

i=i 

is also in involution with the G/s, thereby forming a complete set {Gj}j = o,...2m of commuting 
integrals. Introducing an additional flow parameter ao corresponding to the Hamiltonian Go, 
and letting 

u(a) := u(a) + Godcio, (1-43) 

where 

a:= (a ,...,a 2 m), (1-44) 
we have the following extended system, which is also Frobenius integrable: 

d & Xj = {xj, uj(a)} and d d yj = {yj, uj(ci)} , j = 0, . . . , 2m. (1-45) 

Moreover, since Go is not explicitly dependent on the parameters cij and Poisson commutes 
with all the {Gvj}j=i,...2m> h is in fact a constant. Since all the quantities u,v,qj,pj — > as 
the cij — > oo, this constant value is just 0. 

2. Classical .R-Matrix and Isomonodromy Formulation 

Hamiltonian systems of the type encountered above can be understood within a general 
Lie algebraic setting (cf. [AHP], [H], [RS], [S]) which we now summarize. Let g be the 
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loop algebra of smooth (real or complex), s/(r)-valued functions X(\), defined on a circle S 1 
centered at the origin in the complex A-plane. Define the Ad-invariant scalar product < , > 
on g by: 

< X, y >= -L / tr(X(A)Y(A)) dA. (2.1) 

Z7TZ J s i 

Interpreting this as the evaluation of a linear form X on elements Y G g allows us to identify 
g with an open, dense subspace of the dual space, which will be denoted g*. Under this 
identification, both the adjoint and coadjoint actions of the corresponding loop group Q are 
given by conjugation: 

Ad 5 : X^gXg- 1 , Ad* : X .— gXg~\ geQ, X e g. (2.2) 

As a linear space, g may be decomposed into a direct sum 

g = g++g_, (2.3) 

where g_|_, g_ are the subalgebras consisting of elements X(X) that admit, respectively, holo- 
morphic extensions to the interior and exterior of S , with the elements X e g_ satisfying 
X(oo) = 0. Under the identification g ~ g*, we have 

g+~g-, g-~g+, (2-4) 

where g± C g* are identified as the annihilators of g T . 
Let 

P± 1 1 ^ ^ (2.5) 
P±:X\ — > X± 

be the projections to these two subspaces relative to the decomposition (2.3) (determined, 
e.g., by splitting the Fourier series on 5 1 into positive and negative parts) and define the 
endomorphism R : g — > g as the difference: 

R:=P + -P_. (2.6) 

Then R is a classical R-matrix [S], in the sense that the bracket [ , ]r defined on g by: 

[X,Y} R :=±[RX,Y] + ±[X,RY] (2.6) 

is skew symmetric and satisfies the Jacobi identity, determining a new Lie algebra structure 
on the same space as g, which we denote g#. The Lie Poisson bracket on g* R ~ g* associated 
to the Lie bracket [ , ]r is: 

{/,<7}|x=< [df,dg) R ,X> (2.7) 
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for smooth functions /, g on g R (with the usual identifications, df\x,dg\x G g* ~ g). The 
corresponding group Q R is identified with the direct product Q + x where Q± are the loop 
groups associated to g±. The adjoint and coadjoint actions are given by: 

Ad R (g) :(X + +X_)^ g+X+g- 1 + g.X.gZ 1 

k<? R {g) : (X+ + X_) .— > (g_X + gZ 1 )+ + (g+X_g- 1 )_, (2.8) 
X± e g±, g± e Q± 

Let a := {cii} 1 j L =1 be a set of n distinct real (or complex) numbers and k := {k,ki}f =1 
a set of n + 1 nonnegative integers. Define the subspace g ajK C g R to consist of elements 
X = X + + X- with X + Gg| a polynomial in A of degree k and X_ e g_ a rational function 
of A, with poles of degree {k\, . . . , k n } at the points {ai, . . . , a n }. Then g ajK C g R is Arf^,- 
invariant, and hence a Poisson subspace. The coefficient Y of the leading term in the degree 
k polynomial X + is Ad* R invariant, and hence constant under all Hamiltonian flows. Let 
1 a ,K '■= ^(g*)|ga, K be the ring of Ad* -invariant polynomials on g* ~ g^, restricted to g , K , 
and define, for X e g a , K , the s/(r)-valued polynomial function: 

L(A) := a(X)X(X), (2.9) 

where 



a(A) := JJ(A- ai ) fci . (2.10) 

i=l 

Then the ring X a K is generated by the coefficients of the characteristic polynomial: 

V(X,z) :=det(L(X) - zl), (2.11) 

and the characteristic equation 

V(X,z) = (2.11) 

defines, after suitable compactification (and possible desingularization), an algebraic curve T, 
called the "spectral curve", which is generically an r-fold branched cover of CP 1 . 

It follows from the Adler-Kostant-Symes (AKS) theorem that: 

1) The ring I ayK is Poisson commutative (with respect to the { , } R Poisson bracket). 

2) For a Hamiltonian G G X a ,K, Hamilton's equations are given by: 

HX 

— = {X,G} = [dG + ,X] = -[dG_,X] (2.12) 
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(The last equality following from the fact that [dG, X] = 0, which is the infinitesimal 
form of Ad*-invariance) . 

This means that the spectral curve V is invariant under the flows generated by the 
elements of T a ,n, since any two coefficients of the spectral polynomial V(X, z) Poisson commute. 
It may also be shown [AHP], [AHH1] that, for sufficiently generic initial conditions and 
coadjoint orbits O in g ajK , the elements of I a ,m restricted to O C g a , K , define completely 
integrable Hamiltonian systems; that is, the number of independent generators in X a /t equals 
one half the dimension of the orbit. Through a standard construction [AvM], [RS], [AHH1], 
[AHH2], the flows they generate are shown to linearize on the Jacobi variety associated to V. 

We now specialize to the case {k = = l}j=i,... n - Then X(X) has the form 




(2.13) 



and the Ad^-action is: 




n 



-1 



(2.14) 



where 





n. 



(2.15) 



Hence, each coadjoint orbit in g ajK consists of elements of the form: 




(2.16) 



where 



N :=C+[ l0 ,Y], 



(2.17) 



with Y and C constant matrices, and the iVVs have fixed Jordan normal form. 



Let 




(2.18) 



where 



y := hrY 2 , 



c := tr(YN ) = tr(FC) 



(2.19) 



and 
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Go:=^~<f M 2 (A)^=tr(y^^) + ^triV 2 
711 Js 1 i=l 

i r n 

G l := — f trX 2 (X) dX = aMYNi) + tr(A r AT i ) + ^ 
Kl '- = ^n{ (A " a ^ trX2 ( A ) rfA = T^ 2 ), 



a z - a,- ' (2.20) 



i = 1, ... 77. 

with the integrals j> A=a . taken around a small circle enclosing only this pole. Then y, c and 
{.?Q}i=i,...n are all Casimir invariants and hence constant on each Ad^ orbit, and {G^}i=o,...n 
are independent elements of the Poisson commuting ring I a ,K- Denoting the time parameters 
for the Hamiltonian flows generated by {Gi}i = o,... n as {£i}i=o,...n> the AKS theorem implies 
that Hamilton's equations are of the form: 

dX 

— = {X,G }={(dG ) + ,X] 

M (2.21) 



where 



{X,G l } = -[(rfG l )-,X], i = l,. ..n 



(dG )+ = Y, (dG i )_ = - v ^-, i = l,...n. (2.22) 
Evaluating residues at {A = ai}i=i v .. n and A = oo, we see that this is equivalent to: 

at,- aj — a,- aii ^— ^ — a,- 

i,j = l,...n. 

If we now reinterpret these as nonautonomous Hamiltonian systems, by identifying the 
constants {o^} with the flow parameters {U = ai}i = i v .. n , system (2.22) becomes (cf. [JMMS], 
[JMU]): 

fH^> 

W =[aiy + Wo + E J ! L, JV<] (2.24) 
oa 7 - aj — dj ocLi a\ — a,- 

i,j = l,...n. 
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(Here to has similarly been renamed do, but it does not appear as a parameter in the system, 
and hence represents autonomous flow.) This system is equivalent to the commutativity of 
the following system of differential operators: 



V X := » -XY-N -± * 

OA ^— ' A — di 

i=l 



d N, 



(2.25) 



OCLi A — CL-i 

V :=^-- Y. 
da 

It follows that the monodromy of the operator V\ is preserved under the deformations 
parametrized by {ao, a±, . . . , a n }. 

To obtain the system (1.26)-(1.29), we take the {a.j}? =1 as real, r = 2, n = 2m, Y = 
(so really, we are in the Poisson subspace with k = 0), and choose 

C=(_° 2 f) (2.26) 

and the iVj's with Jordan normal form ^ . The Ad^ orbit O C g ajK then consists of 
elements of the form: 

( -XiVi -yl \ 



2 / 2^ A-a 

7 2=1 

Here, we have parametrized the orbit by elements of Mq := R 4m /(Z2) 4m , where 

((xi, . . . , i2 m ), (yi, . . . , y2m)) € R 4m , and the (Z 2 ) 4m action is generated by reflections in the 

coordinate hyperplanes: 

(ei,...,e 2m ) : ((xi, . . . , x 2m ), (yi, • • • , y2m)) 1 — > ((ei^i, . . . , e 2m a;2 m ), (eiyi, . . . , e 2m y2m)) 

{fi = ±l}i=l,...2m- 

(2.28) 



The map J : M — > g* R defined by 



. XiVi -Vi 

2m I 9 



-2 0~) + 2^ X-a % 1 ^ 2 ' 29 ^ 

is a symplectic embedding with respect to the standard symplectic structure: 

2m 

uj = ^2dxi Adyi (2.30) 
i=i 
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on M and the orbital (Lie-Kostant-Kirillov) symplectic structure on Oq. (This is a special case 
of the moment map embeddings into rational Poisson subspaces of the dual of loop algebras 
developed in [AHP], [AHH3].) Evaluating the Hamiltonians {Gi}i=i,...2m defined in (2.20) 
for this case gives those of eqs. (1.26) and the {%}, j ^ components of the equations of 
motion (2.24) are equivalent to (1.29), up to quotienting by the (Z 2 ) 4m action (2.28). 

To obtain the system (1.40)-(1.45), we again take r = 2, n = 2m, but choose 
with the Ni's again of Jordan normal form ( ^ J. In this case, the Ad^. orbit 0\ C g ajK 



consists of elements of the form: 



. -XiVi -yf 



N ' i=l 

Here, we have parametrized the orbit by elements of M\ := R 2 x R 4m / (Z2) 4m , where (xo, yo) £ 
R 2 , ((xi, . . . , X2m), (yij • • • , y2m)) £ R 4m , and the (Z 2 ) 4m action is again given by (2.28). The 
map J\ : Mi — ► g|j defined by 

( -XiVi -vl 



( v —- + — \ 1 

Ji : ((x , xi, . . . , x 2m ), (y , yi, • • • , y2m) 1 — M _ 2 2 _ yo 2 J + 2 ? 

\ / A — 1 



i=l A - a * 



(2.33) 

is again a symplectic embedding, now with respect to the symplectic structure: 



2m 



to = dxo A dyo + A (2.34) 



i=l 



on Mi and the orbital (Lie-Kostant-Kirillov) symplectic structure on 0\. Evaluating the 
Hamiltonians {Gi}i=o,...2m defined in (2.20) for this case gives those of eqs. (1.41), (1.42) 
and the equations of motion (2.24) are equivalent to (1.45) - again, up to quotienting by the 
(Z 2 ) 4m action (2.28). 

For this case, eqs. (2.19), (2.31), (2.33) give y = 0, c = 1 and tr(iV 2 ) = 0, so 

-det(X) = itrPO = A + G„ + £ ^- := |W (2.35) 
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where P(X) is a monic polynomial of degree 2m + 1 



P(\) = A 2m+1 + P 2m A 2m + ... + P 



(2.36) 



with 



2m 2m 2m 



Pim = G - } a,i, Pim-\ = > Gi - G } at, etc. 



(2.37) 




From eqs. (2.9)-(2.11), we have the following equation for the spectral curve F: 



V{\z) = z 2 - d{X)P{\) = 0, 



(2.38) 



which shows that it is hyperelliptic, of genus g = 2m, with 2m + 1 of its branch points located 
at {ai, . . . dim, 00 } an d the others determined by the values of the roots of P(X). For the 
system (1.26)-(1.29) associated to the sine kernel, the leading term in (2.18) vanishes, so the 
polynomial P(X) is of degree 2m and the curve has genus g = 2m — 1, with no branch points 
at oo. For the autonomous system (2.23), these curves are invariant under the flows. For the 
nonautonomous system (2.24), which in these cases reduces to (1.29) or (1.45), it would be 
of interest to determine the dependence of the spectral curves on the deformation parameters 
{cii}i=o,...2m through a system of PDE's involving the spectral invariants {Gi}i = o,...2m alone. 
For the single finite interval case, this is given in [JMMS], [TW1] by the r-function form 
of the Painleve equation Py, while for the single semi-infinite interval case, this is given in 
[TW2] by the corresponding r-function form of Pjj. 
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